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It is well known that the concept of return difference plays a central role 
in the classical theory of linear feedback systems developed by Black, Ny- 
quist, Blackmail and Bode. This concept, which relates to the influence of a 
single, algebraic system-constraint of the form /•> = 7/1 where f\ and f 2 re- 
spectively may be regarded as a controlling signal and a controlled signal, 
retains its prominence in the subsequent signal-flow graph theoretic exten- 
sions by Mason. It is particularly pertinent to the study of the stability of 
the system, its degree of immunity from parameter variations, and the de- 
termination of its transmission ami driving-point properties. 

This paper reports on a generalization of Blackman's equation and on 
some generalizations of Bode's return difference theorems. Here attention is 
focused not on a single constraint of the form f 2 = 7/1 , but instead on a set 
of constraining equations. 

I. INTRODUCTION 

It is well known that the concept of return difference plays a central 
role in the classical theory of linear feedback systems developed by 
Black,' Nyquist, 2 Blackmail 3 and Bode. 4 This concept, which relates to 
the influence of a single algebraic system-constraint of the form/ 2 = 7/1 
where f\ and / 2 respectively may be regarded as a controlling signal and a 
controlled signal, retains its prominence in the subsequent signal-flow 
graph theoretic extensions by Mason. 5 It is particularly pertinent to the 
study of the stability of the system, its degree of immunity from parame- 
ter variations, and the determination of its transmission and driving- 
point properties. 

This paper reports on a generalization of Blackman's equation 3 and 
on some generalizations of Bode's return difference theorems. 4 Here at- 
tention is focused not on a single constraint of the form f 2 = 7/1 , but 
instead on a set of constraining equations. For this reason the results 
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are applicable to systems which possess in the usual physical sense a 
multiplicity of feedback loops. However it will become clear that the 
results are not restricted to situations of this type. 

In Section II we describe the basic system considered throughout the 
paper. Section III presents an explicit expression for w(X), the trans- 
mission or driving-point function to be studied. We then introduce defini- 
tions, based on a simple topological characterization of the relation be- 
tween the system input and output variables, of the loop-difference 
matrix, the null loop-difference matrix, and the complementary loop- 
difference matrix. The determinants of these matrices are of fundamental 
importance in the subsequent discussion. In Section VI we derive a gen- 
eralization of Blackmail's interesting equation. The material in Sections 
VII and VIII relates to generalizations of Bode's well-known return 
difference theorems. Some applications and illustrations of the theory 
are discussed in Section IX. 

II. THE BASIC SYSTEM AND THE SET OF EQUATIONS g 

We shall be concerned throughout with the transmission and driving- 
point properties of the structure shown in Fig. 1, an arbitrary linear 
time-invariant two-port network containing no independent sources. 

The Laplace-transformed equilibrium equations which implicitly de- 
fine the external properties of the two-port are of the form: 

E («*** + Mi) =0 (ft = 1, 2, • • • , K) (1) 

j 

where ej and ij respectively are branch voltages and currents and the 
cthj and fikj are functions of the complex-frequency variable. 

We wish to focus attention on the influence of a prescribed subset of 
the linear constraints implied by (1). This subset is assumed to be ex- 
pressible as 

/*= E ykifi (ft = 1,2, ••-,?) (2) 

in which each fi is one of the e ; or one of the ij . It is convenient to 
interpret these relations as corresponding to a set of q controlled (i.e., 




Fig. 1 — Two-port network. 
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dependent) sources, where the quantity on the left-hand side of each 
equation is a controlled-source output and each /,• on the right-hand 
side is a controlling variable. Thus we take the viewpoint that the two- 
port network contains a controlled-source subnetwork characterized by 
(2) even when the two-port does not contain devices such as vacuum 
tubes or transistors with which controlled-sources are ordinarily asso- 
ciated. We shall denote by g the set of controlled-source constraints (2). 
Let the components of a g-vector *F be the f k (k = 1, 2, • • • , q) ar- 
ranged in any one of the ql possible orders, and let the components of a 
p-vector 4> be the controlling variables/, (i = q -f 1, q + 2, • • • , q + p) 
arranged in any one of the p! possible orders. Then (2) can be written 
as W = X4>. We shall call the q X p matrix X the controlled-source 
matrix or the matrix of controlled-source coefficients. 

III. EVALUATION OF THE TRANSFER AND DRIVING-POINT FUNCTIONS FOR 
THE TWO-PORT 

Our primary concern here is the determination of the influence of g 
(i.e., of the controlled-source subnetwork) on the transfer and driving- 
point functions for the two-port network. 

Let ?/ 2 = w(X)yi represent the relation to be studied in which y 2 
and 7/1 respectively are response and excitation functions (each function 
may be either a voltage or a current). If w(X) is a driving-point immit- 
tance, y 2 and //i respectively are a current and a voltage or a voltage and 
a current at the same port. 

Consider now an evaluation of w(X). With ?/i and the components 
of W treated as independent variables, we apply the superposition 
theorem to obtain f 

2/2 = d !h + BW (3) 

where d and B are defined by the equation. In particular when X = 0, 
W = 0, and ?/o = diji . That is, d = w(0). 
Similarly we can express <P as 

* - A//, + OF (4) 

where the matrices A and C are defined by the equation. After pre- 
multiplying both sides of (4) by X and using W = X<I>, we find that 

W = [1, - XCF'XA//! (5) 

where \ q is the identity matrix of order q. 

t We are assuming that 1/2 is uniquely determined by ?/i and*. See the relevant 
discussion in the next section. 
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Fig. 2 — Basic flow-graph representation of the relation between v/i and ;/■•. 

From (3), (5), and d = w(0) 

w(X) = ytiyxT 1 = w(0) + B[lg - XCf'XA. (6) 

IV. THE BASIC FLOW GRAPH 

In order to express w(X) in terms of quantities that are related to 
concepts of importance in the classical theory of linear feedback systems, 
we first represent the relation between yi and y 2 by the signal-flow 
graph 5 in Fig. 2. The matrix equation associated with "node-group" 
«I> is (4), and so forth. This topological characterization of the influence 
of 5 plays a central role in the subsequent discussion. We shall call 
A, B, C, and d "flow-matrices." 

It is worth emphasizing at the outset that although the signal-flow 
graph exhibits feedback in the sense that the outputs of the controlled 
sources (i.e., the elements of W) influence the value of the controlling 
variables (i.e., the elements of <I»), it is not necessary that feedback 
exist in the network in the usual physical sense. For example, the 
physical system may comprise simply a set of driving-point impedances 
and the controlled sources may represent a certain subset of these 
one-port elements. Feedback arises in the topological characterization 
merely because of the form of the equations that we have chosen to 
write. 

Before proceeding it is important to note that for some choices of g 
it may not be possible to characterize the two-port by equations of the 
form (3) and (4) and hence by the flow graph in Fig. 2. The superpo- 
sition theorem implies that when ?/•> and <I> are uniquely determined by 
yi and W, the relations are of the form (3) and (4); it does not imply 
that such relations exist. The one-port network shown in Fig. 3 is one 
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Fig. 3 — Single controlled-source one-port. 

of the simplest structures which illustrate the difficulty. For this net- 
work, with w(X) the well-defined driving-point admittance, all four 
flow matrices A, B, C, and d fail to exist. 

There are much more sophisticated situations of this general type. 
Nevertheless, almost all feedback networks of interest do not exhibit 
this degeneracy. Furthermore when the difficulty does occur, it is 
generally possible to consider the limiting form of a network for which 
the degeneracy is not present. For example, if the one-port in Fig. 3 
is modified by adding a series resistor, the flow matrices A, B, C, and d 
become well-defined and the expression obtained for the driving-point 
admittance reduces to (X) -1 as the value of the series resistor approaches 
zero. 

With this motivation we state the 

Assumption: It is assumed throughout that the relation between iji and y 2 
can be represented by a signal-flow graph of the type shown in Fig. 2. 

Consider now the definition of three matrices that are related to the 
flow graph. 

4. l The Loop-Difference Matrix for the Branch X: Fj(Xi) 

The loop-difference matrix for the branch X is defined as follows. 
Introduce an additional node-group P in the flow-graph of Fig. 2 by 
replacing X with a cascade of any two branches Xi and X 2 such that 
X,X 2 = X. Let the orders of X! and X 2 respectively be q X m and m X p. 
Next, break the feedback loop by splitting P into a source node-group 
P' and a sink node-group P" to obtain the graph in Fig. 4. 

With yi set equal to zero, suppose that an arbitrary signal vector 
S H1 , of order m, is applied to P\ The resulting signals at node-groups 
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Fig. 4 — Signal-flow graph for defining the loop-difference matrix and the null 
loop-difference matrix. 

W, <P, and P" respectively are XiS„, , CXjS,,, , and XoCXiS™ . The vector 
signal difference at P',P" (i.e., the signal at P' less that at P") is [1„, — 
X 2 CXi]S m . We shall call [1,„ — X 2 CXi] the loop-difference matrix for 
the branch X and denote it by F s (Xi), which indicates explicitly its 
dependence on X t . When Xi , X 2 , and C are scalars, Fg(Xi) becomes 
independent of Xi and reduces to Mason's flow-graph definition of the 
loop-difference for the branch X. 

4.2 The Null Loop-Difference Matrix for the Branch X: Fg(Xi) 

This matrix is evaluated under the condition that y± is adjusted so 
that 2/2 = 0. Specifically, if an arbitrary signal vector S m is applied to 
P' in Fig. 4, the signal reaching y-> by way of the branch B is BXjS,,, 
so that, if 7/2 is to be zero, //i must be — </ _I BXiS„, . The total signal 
arriving at 4» is — d~ l ABXiS m + CXjS m , and hence the signal at P" 
is [X 2 CX! - rf" 1 X 2 ABX 1 ]S,„ . Thus the vector signal difference at P',P", 
under the condition that yi is adjusted so that ?/ 2 = 0, is [l m — X 2 CX! -+- 
d" 1 X 2 ABX 1 ]S m . We call [l m - XoCXi + (/"'XoABXj] the null loop- 
difference matrix for the branch X and denote it by F 5 (Xi). When A, 
B, C, Xi , and X 2 are all scalars, Fg(Xi) reduces to Truxal's definition" 
of the null loop-difference for the branch X, and is independent of the 
choice of Xi . It is convenient to write: Fg(Xi) = [1„, — X 2 CXi], where 
C = C - f/ _1 AB. 

4.3 The Complementary Loop-Difference Matrix for the Branch X: Fg(Xi) 

The complementary loop-difference matrix, denoted by Fg(Xi), is 
defined by the requirement that when an arbitrary S m is applied to 
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P', F 8 (X])S m is the vector signal difference at P',P" under the condition 
that yi is adjusted so that (?/i + y 2 ) = 0. This matrix is therefore ob- 
tained from the expression for the null loop-difference matrix by setting 
d = 1. That is, FbCXO = [l,„ - X 2 CXJ where C = C - AB. The comple- 
mentary loop-difference matrix is of utility when d = 0. In such cases 
the null loop-difference matrix is not defined. 

It is evident that the concepts of loop-difference matrix, null loop- 
difference matrix, and complementary loop-difference matrix need not 
be restricted to the branch X in Fig. 2; they relate without ambiguity 
to any branch in a signal-flow graph which possesses a single input node 
and a single output node. 

4.4 Circuit Interpretations of C, C, and C 

It is important to note that the matrices C, C, and C possess very 
explicit circuit interpretations. When the elements of W (i.e., the outputs 
of the controlled sources) are treated as independent variables ,_ the p- 
vector of voltages and currents * is equal to CW, CW, and C*F re- 
spectively when yi is set equal to zero, yi is adjusted so that y 2 = 0, and 
2/1 is adjusted so that (3/1 + 1/2) =0. The evaluation of these matrices 
from the circuit is simplified considerably by the fact that the controlled 
sources are treated as independent sources. 

v. the determinants det F a (Xi), det £ s (Xi), and det Fg(Xi) 

Our primary interest in the matrices F 8 (Xi), Fg(Xi), and Fg(Xi) is 
with regard to their determinants. Here we wish to establish two ele- 
mentary properties of these determinants which add to an understanding 
of the character of the expression for w(X) presented in the next section. 
Further properties of these determinants are considered subsequent to 
the derivation of the expression for w(X). 

Although the matrices F«(Xi) and F 8 (X0 generally depend upon 
the choice of Xi , it is true that 

Lemma I: 

det F*(Xi) = det Fx(l g ) 

dett s (Xi) = dett 9 (l 9 ). 

The lemma is proved in Appendix A. It implies, of course, that the 
determinants are independent of the_choice of X! . This property is 
evidently shared by det F B (Xi) since P B (Xi) can be obtained formally 
from $»(Xi) by setting d = I. 
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Lemma I is often of assistance in evaluating the determinants. 
For example, if the normal rank of X is unity, as is the case when $ 
contains a single equation, X can be written as XiX 2 where X x and X-. 
respectively are q X 1 and 1 X J) matrices. In such instances the de- 
terminants of the loop-difference matrix, null loop-difference matrix, 
and complementary loop-difference matrix can be expressed respectively 
as simply 1 — XoCXi , 1 — X 2 CXi, and 1 — X 2 CXi . Similar simplifica- 
tions can of course be exploited if C, C, or C is of unit normal rank. 

With the result stated in Lemma T as motivation, we shall throughout 
the remainder of the paper denote det Fg(Xi), (let F^XO, and det 
Fg(Xi) respectively by det F 8 , det F 8 , and det F 8 . 

Recall that the components of W and <I> respectively are the con- 
trolled variables and controlling variables arranged in any definite 
orders. Thus the matrices X, F s (Xj), and F R (X,) are not uniquely 
defined by the details of the two-port structure and the set of controlled- 
source constraints 5- Nevertheless, 

Lemma II: det F R and det F 3 are invariant with respect to the ordering of 
the components of W and <P. 

In other words, the determinants are uniquely defined by the details 
of the two-port and $• The proof of this important result is straight- 
forward: Let <!>' and W' respectively denote vectors obtained from <I> and 
W by any reordering of the components. Then <I>' = U* and W' = V*F, 
where U and V are nonsingular matrices. In terms of «I>' and W' the 
controlled-source constraints read W' = X'*' where X' = VXU~ , and 
the equations corresponding to (8) and (4) are 

//, = dyi + B'n" = d !h + BV^r' (7) 

«>' = a'//! + C'*r' = UA//, + UCrV. (8) 

Finally, note that for all d ^ 

det [1 3 - X'(C - d-'A'B')\ = det [1„ - VXU VUCV -1 - rT'UABV 1 )] 

= det [1, - X(C - f/ _1 AB)]. 

VI. GENERALIZATION OF BLACKMAN's EQUATION' 

At this point we are in a position to state and prove the following 
generalization of Blackman's classical result. 



Theorem I: If w(0) ^ 0, 



•«—<«£$■ 
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6.1 Proof 

Consider equation (6), which for convenience is repeated below: 

w(X) = VtM' 1 = w(0) + B[l a - XCf'XA. (9) 

Recall that adj (Iff), the adjoint of an arbitrary square matrix M of 
order q, is defined by 

M adj flff) = adj (M.)M = l q det M. 

Thus w(X) can be expressed as 

w(X) = w(0) + B adj (l a - XC)XA[det (1, - XC)]" 1 

" det [1, - XC] + wjoy'B adj (1„ - XQXA l 
det [1, - XC] J ' 



= w(0) 



Since det Fg = det [l fl - XC], we must prove that det [1, - XC] + 
m)(0) _1 B adj (1„ - XC)XA = det F 5 or, more explicitly, that 

det [l f - XC] + w(0)- l B adj (1, - XOXA 

= det[l g - X(C - w(0) -1 AB)]. 

We first prove the following result which will be used also in a later 
section. 

Lemma III: Let G and H respectively be non-singular matrices of orders 
n and m, and let I and J respectively denote matrices of orders n X m 
and m X n. Then, 

det[G + IH^Jjdct H = det\R + JG -, I]de/ G. 

The proof follows from the lemma established in Appendix A which 
states that det [1„ + DEJ = det [l m + EDI for arbitrary matrices D 
and E respectively of orders n X m and m X n. Taking D = G 'I and 
E = H _1 J, we have: det |1„ + G 'lH"'j] det H det G = det [l m + 
H'jG'l] det G det H. Moreover, det [1„ + G~'IH _1 J] det G = det 
[G + IIT'J], and det [1„, + H 'JG ! I] det H = det [H + JG -1 I]. 

The identity (10) is a direct consequence of the following corollary 
of Lemma III. 

Lemma IV: Let G, I, and J respectively denote matrices of orders n X n, 
n X 1 , and 1 X n. Then, 

del[G + IJ] = det G + J adj(G)l. 

To establish this result, consider Lemma III with H = 1 and m = 1. 
Clearly when G is nonsingular, det [G + IJ] = (1 + JG _I I) det G = 
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Fig. 5 — (a) Network with the open-circuit transfer impedance from port 1 
to port 2 equal to w(X); (b) modified network with driving-point impedance 
w(X). 

det G + J adj (G)I. A continuity argument of the type used in Ap- 
pendix A shows that Lemma IV is valid also when G is singular. 

Expression (10) can be obtained from Lemma IV by taking 
G = (1, - XC), I = XA, and J = u>(0) -1 B. This completes the proof 
of Theorem I. 



6.2 Remarks Relating to Theorem I 

When w(X) is a driving-point impedance function, the null loop- 
difference matrix is evaluated under the condition that the current at 
the terminal pair at winch the impedance is denned is adjusted so that 
the voltage at the terminal pair is zero. Thus, in this case, the null 
loop-difference matrix is equal to the loop-difference matrix evaluated 
from the network when the terminal pair is shorted. Similarly when 
w(X) is a driving-point admittance function, the null loop-difference 
matrix is equal to the loop-difference matrix evaluated from the network 
when the terminal pair is open-circuited. 

The conceptual and computational simplifications mentioned in the 
last paragraph are not valid if iv(X) is a transfer function. f However, 
by introducing an additional controlled source, the evaluation of the 
open-circuit transfer impedance or short-circuit transfer admittance 
of any three-terminal network can be reduced to the evaluation of a 
driving-point function. This reduction is illustrated in Fig. 5 for the 
transfer impedance case. 



t Of course we assume here that w(X) is not both a transfer function and a 
driving-point function associated with the two-port in Fig. 1. 
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6.3 Relation of Theorem I to Earlier Work 

The expression given in Theorem I reduces to Blackman's classical 
result 3 when it is applied to the driving-point immittance case in which 
g characterizes a single unilateral amplifier, and the equation is expressed 
in terms of physically-defined open-circuit and short-circuit return dif- 
ferences. 

A particularly succinct derivation of an equation of this form (for 
the single amplifier case) was subsequently presented by Bode 4 who 
exploited the properties of the determinant of the network immittance 
matrix. In addition Bode introduced the useful concept of return differ- 
ence with respect to a two-terminal element. Signal-flow graph interpre- 
tations of Blackman's expression were later considered by Mason 5 and 
by Truxal. 6 Our proof of Theorem I is a generalization of Truxal's work. 

Exploiting a suggestion by Bode, 4 Mulligan 7 has stated an entirely 
different generalization f of Blackman's equation for cases in which a 
multiplicity of nonreciprocal elements are to be considered. His result is 
an explicit expression obtained by repeated application of Blackman's 
original result. 

G.4 Generalization of Bode's Relation between Feedback and Impedance 

Theorem II: Let del F' s denote the determinant of the loop-difference matrix 
for the set of controlled-source constraints g, evaluated under the condition 
that Hi = 7/' _1 //2 . Let u»(X) — «'(0) not vanish identically in X. Then 
del F'g is a linear-fractional function of 77 that vanishes in t\ if and only if 

1) = W>(X). 

The proof is based on Lemma IV: From the signal-flow graph obtained 
from Fig. 2 by adding a branch from y 2 to y x with transmission if\ 
it is clear that 

detF' 5 = det {1, - X[C + AB(t, - d)" 1 ]}. 
Using Lemma JV, 

detF' 5 = det [1, - XC] - B adj (l g - XC)XA(ij - rf)" 1 

Thus the equation det F' s = implies that r\ = w(X), assuming that 

f Subsequent to preparing I his paper and soliciting comments from colleagues, 
it came to the writer's attention thai a result similar to Theorem I is contained 
in the recently published book by Yutzc Chow and Etienne Cassignol: Linear 
Signal -Flow firaphs and Applications, John Wiley and Sons, New York, 1962. 
They consider the situation that corresponds here to the special case in which J? 
is a set of equations of the form/* = ti/a-h? (k = L 2, . . . , q). Their proof is con- 
siderably different from ours. 
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B adj (lg — XC)XA = w(X) — w(0) does not vanish identically in X. 
Of course when w(X) is a driving-point impedance function, the 
condition i/i = tT l y% corresponds to adding an impedance — r\ in series 
with the one-port. 

6.5 The Expression for w(X) when w(0) = 

Theorem I does not apply when w(0) = 0. In fact, in all such instances 
the null loop-difference matrix does not exist. The following corollary 
which involves the complementary loop-difference matrix is of assistance 
in these cases. 



Corollary I: If w(0) = 0, 

w(X) = 



del F g - del F 8 
detFx 



This result is readily established: Here d = in the flow graph of Fig. 2. 
Suppose that the branch d is replaced with one having unit transmission. 
Then the ratio of ?/ 2 to ?/i for the resulting graph is w(X) -\- 1. However, 
using Theorem I and the definition of Fg, w(X) + 1 = det Fg[det F 8 ] _1 . 

VII. THE MATRICES F 58o AND F 55o , AND THE NOTION OF A "RESIDUAL SET 
OF EQUATIONS" OBTAINED FROM ft 

In this section we generalize the definitions of the loop-difference 
matrix and the null loop-difference matrix and then introduce the notion 
of a "residual set of equations" obtained from $, in order to both 
facilitate and provide a firm general basis for the subsequent derivation 
of some fundamental properties of det F 8 and det F 8 . The complemen- 
tary loop-difference matrix need not be considered separately since 
it is merely a special case of the null loop-difference matrix. The ma- 
terial presented in this and the next section is in many respects a gen- 
eralization of Bode's classical theory relating to return difference with 
respect to a single element. 

7.1 The Matrices F 55o and F 5Bo 

The relation between y x and the quantities O, W, and y-> clearly 
remains unchanged if the signal-flow graph shown in Fig. 2 is replaced 
with the graph in Fig. 6, where X is an arbitrary q X p matrix. Let 
the loop-difference matrix and null loop-difference matrix for the branch 
(X — Xo) in Fig. 6 be defined in the same manner as for X in Fig. 2, 
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Fig. (\ — Signal-flow graph for the definition of F 33o and F 33o . 

and let $o denote the "reference set" of equations obtained from (2) 
when the controlled-source matrix X is set equal to X . We introduce 

Definition I: The matrices F 33o and F 33o respectively denote the loop- 
difference matrix and the null loop-difference matrix for the branch (X — X ) 
in Fig. (5. When X = 0, we write: F 33o = F 3 , F 33o = F 3 . The determi- 
nants del F 33o and del F 33o respectively are referred to as the determinant 
of the loop-difference matrix, and the determinant of the null loop-difference 
matrix for the set of controlled-source constraints %, with respect to the 
reference set jyo . 

The branch (X — X„) in Fig. 6, which corresponds to the matrix 
relation M*' = (X — X )4», may be interpreted as characterizing a set 
of controlled sources which focuses attention on the departure of the 
elements of X from those of X or, equivalently, on the departure of % 
from Jo . 

We wish to prove 

Theorem III: Let del F 3 | 3=3l) and det F 5 | 3=3o respectively denote the 
determinants obtained from det F 3 and det F 3 by replacing each element of 
X with the corresponding element of X . Then, 

1tV det F* 

det F 33o = 



det F 33o = 



det F 3 | 3=3o 

det F 3 

det F 3 | 3=3o 



Node-group *F in Fig. can be eliminated to obtain the flow graph in 
Fig. 7. For example, the branch transmission from W' to 4> in Fig. 7 
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Fig. 7 — Flow-graph obtained from Fig. 6 by eliminating node-group W. 
is the transmission of the subgraph in Fig. 8. From Fig. 7 it is clear that 



Thus, 



detF 3Bo = det [1, - (l p - CX ) _1 C(X - X n )]. 



detF 33o = det [(1, - CX ) _1 ] det [1, - CX] 



det F 3 
det F s | 3=3o 

Now consider det F 33n . Since 



w 



(X) - W (0) **** = w (Xo) ^^ 

det F 5 det F 53o 



= [d + BX (l„ - CXn)-^] 



det {W det F* |«_. 



35o 



8 IS— So 



fll) 



(12) 



(13) 



(14) 



det F 3 

where (12) has been used to obtain the last expression in (14), we have 

det F B = [1 + d _1 BXo(l p - CXo) _1 A] det F 33o det F 3 | 3=3o (15) 

- {det [l p - CXo] + rf _1 BX adj (1, - CX )A| det F 33o . (16) 

Although it is assumed in (13) that det F 3 and det F 33o do not vanish 




!q C 

Fig. 8 — Subgraph for evaluation of transmission from V' to «i> in Fig. 7. 
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identically in the complex-frequency variable, it can readily be shown 
with a continuity argument that the validity of (16) does not depend 
upon these assumptions. According to Lemma IV, the expression within 
the branches in (16) is equal to det [l p — CXo]. Therefore 

det F* 



det F 53o = 



det F s | 5=3o 



7.2 The "residual set of equations" 

Recall that the set g is the collection of q equations: 

A=EW, (Jb-l,2,...,ff). (17) 

»-«+i 

Let (R denote an arbitrary subset of (Ro , the set of all ordered pairs of 

integers (j, i) such that 1 ^ j ^ q, q + 1 ^ i ^ q + p. Let E denote 

in 

a sum over all integers i such that (k,i)e(R [i.e., such that (k,i) is an 
element of (Ft], and denote by (R* the complement of (R with respect to 
(Ro . Then it is certainly true that 

A - /*' + E yuft (* = 1, 2, • • • , q) 

61* 

where 

ti = EW, (* = 1,2, •■■ ,q). (18) 

(R 

In accordance with the controlled-source interpretation of (17), the 
f k ' are the contributions to the controlled-source outputs associated 
with the subset of coefficients J7*, | (fc,i)e(R}. Let %-Gl denote the set of 
/• equations (r ^ q) obtained from (18) by omitting all equations of the 
form // — 0. The set % • (R is referred to as a residual set of equations 
obtained from %. 

VIII. THEOREMS CONCERNING det Fg.® AND det Fg.fl 

To each choice of (R there corresponds a signal-flow graph characteri- 
zation of the relation between \j\. and \ji of the type shown in Fig. 2 
where the elements of the flow matrices analogous to A, B, C, and <7 
are independent of those literal coefficients y k i for which (fc,i)e(R, and 
where the branch analogous to X is associated with the set of equations 
iy-(ft. Accordingly, each choice of (R defines a set of controlled-source 
interpretablc equations $-<R, a pair of determinants det Fg.^ and det 
Fg.(R , and an initial state of the physical system [i.e., a system obtained 
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by setting 7*, = for all (k,i) contained in a]. The primary purpose 
of this section is to derive some fundamental results that relate det 
F 5 . (ltl and det fVoti respectively to det Fa-,,,, and det Fg.m, , where (Ri 
and (R 2 are any two subsets of fll . The set g-(Ro = $ is of interest 
here only in that gflh and 5-0*2 are obtained from it in a prescribed 
manner. 

Our first objective is to relate det F^ and det F 5 . -t respectively to 
det F 5 and det F 5 . Consider 

Lemma V: Let A', B', C, and cV respectively denote the flow matrices in 
Fig. 7 which correspond to A, B, C, and d in Fig. 2. Let *, W', A', B', 

and C be partitioned as follows: 



<I> 



[3 



(p - .)' 

B' = [B/ B 2 '], 

r (q — r) 



\jr' 



m 



\r 

(q - r] 



A' = 



U'J( 



p - s) 



C = 



sY 



C n C 12 ]s 

r (q — r) 

and let all of the nonzero elements of (X — X ) be restricted to X rs , the 
r X s subm-atrix standing in the upper left-hand corner. Then det F BBo 
and det F 3 g respectively are equal to the determinant of the loop-difference 
matrix and the determinant of the null loop-difference matrix for the branch 
X r „ in the flow graph of Fig. 9. 

The proof follows at once from the expressions for det F SS() and det 
F 85o in terms of d' and the submatrices of A', B' and C. The details 
are omitted. 

Suppose now that the elements of W and 4> are chosen so that all of 
the coefficients ta-JCA'^cGI] are contained in the r X s submatrix in the 




Fig. 9 — Flow-graph relevant to Lemma V. 
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upper left-hand corner of X, where s is the smallest integer for which 
this is possible, and that X Q in Fig. 7 is obtained from X by replacing 
the 7*»[(fc,i)e(R] by zeros. It then follows from Lemmas II and V that 
det F 35o = det F 8 .<r and det F 35o = detF 5 . m . Thus a direct application 
of Theorem III at this point proves 



Theorem IV. 



det F*, m = 



det F a 



det Fg 



yti = 

(k,i)m 



rfc/Fg.a, - 



det F B 



det Fa 



7k = 

(k,i)e6l 



Note that det Fg.m and det Fg.^ are linear-fractional functions in each 
of the y ki [(k,i)e6{*]. This fact is often of assistance in evaluating the 
determinants. 



Now consider 



det Fa 



y ki = , det F 3 

,,.»„, and 
(k,i)e(R 



y k i = 
(&,t)e<R. 



Recall that these quantities are respectively the determinant of the loop- 
difference matrix and the determinant of the null loop-difference matrix 
for the branch X in Fig. 2 when y ki = 0[(k,i)e(R]. Consequently it 
follows from Lemma IT, Lemma V, and the significance of %•<&* that: 



det F 5 



y ki - _ detFg.^. 
(k,i)e6l 



dot Fg 



7*. = 
(k,i)eQ\ 



where 



det Fg., 



7ki = 
(fc,t)e(R 



and 



det Fg. m . 
det F B . m . 
det F 5 . ( „. 

det Fg. m 



yki = 
(k,i)etil 

y k! = 
(k,i)e<R. 

yki = 
(fc,«)e(R 

| (fc,?')e(R. 

7" = 
(fc,t)e(R 



(19) 



(20) 



respectively 



are equal to det Fg.^. and det Fg.(R. evaluated from the flow graph or 
directly from the circuit under the condition that 7*; = for all (k,i) 
contained in (R. Theorem IV and identities (19) and (20) imply 
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Theorem V: Let (Ri and (R 2 denote two arbitrary subsets of (Ro . Then 



det Fg.gti _ 

r/e£F s .(Ri 



7ki = 



7k = 
(fc,t)e<Ri 

7* = 



detPfrfa det f *.<&,• 



7m = 

(fc,i)e(Ri 



We wish now to focus attention on the particular situation in which 
(Ri and (R2 are disjoint. We shall prove 

Theorem VI: Let (Ri and (R 2 be any two disjoint subsets of (R . Then 



det Fg.R, 
de<F B .(R 2 de«F 8 .(R 2 

deiFjj.Kj 

^ Fg-flt, = 

dei F*. 



7M = 
(k,i)e(Ri 



3-(R 2 



dei F«. 



S-012 



7*.- = 
(fc,i')e(Ri 

7*,- = 
(/c,?')e(R 2 



7a- = 
(k,i)e(Ri 



Consider Theorem V. Observe that here g-(Ri and g-(R 2 respectively 
can be regarded as residual sets of equations obtained from ^ ■ ^2* and 
5*(Ri*. Formally, in accordance with the notation introduced earlier, 
% • (Ri = ( g • (R 2 *) • (Ri and $ • (R 2 = ( g • (Ri*) ■ (R 2 . Using Theorem IV 



det Fg.ffl, = 



detF 



8-(R 2 * 



det F 5 . m2 



7M = 

(fc,z')e(Ri 



Hence 



detF: 



S.(R2* 



7m = 

(k,i)e(R 2 



det Fg.^j 



7m = 

(A;,Oe(Ri U (R 2 



= det F 



S.(Ri 



= det F, 



3.(51] 



7*< = 

(ft,i)e(R2 

YM = 

(/e,i')e(R 2 " 



(21) 
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Next let (R 3 = (<JU U (R 2 )* and note that g-(R 3 = (g-Ou*)-^ 
= ( 5 • <R2*) • (Rs • Again using Theorem IV 

det Fg.ffl,. _ det F 8 .<r,. 



detF 



B(Rj 



det Fb.^. 



7w = 

(fc,i)e(Rs 



det F b .(r 8 . 



7ki = 
(A;,t)e(5l3 



Thus 



det F g . ( 



y ki = = det Fg-flt, 

(Jb,t)e<Ri U (R 2 



T*< = 
(fc,i)e(Ri U (R 2 



= detF 3 . 



Slj' 



(22) 



7*i = 

(k,i)e(Ri U (R 2 . 

Finally, from (21), (22), (19) and the fact that (21) remains valid 
when the subscripts 1 and 2 are interchanged, it follows directly that 



det F 



8.<Ra" 



detF 



5. (Hi 



in = 
(k,i)e6\ 2 



det Fg.at! 



det F 



3.m 2 



TW = 
(fc,i)e(R 2 



7« = 
(fc,i)e(Ri 



7*< = 
(fc,i)e(Ri 

when (Rj and (R 2 are disjoint. It is obvious at this point that a similar 
argument suffices to establish the second identity stated in the theorem. 
It is worth stating explicitly the following direct corollary of Theorem 
VI. 

Corollary II: Let N\ and No be two arbitrary disjoint subsets of {1,2, • • • , q\. 
Let $1 and fjj respectively denote the subsets of equations obtained from % 
by including only those for which keNi and keNo . Then 



del F gl 
det F 52 



det F* 



det F 5l 



det F g , 
det F 8l 



T« = 

keN 2 



det Fg 2 det Fg 2 



7w = 

keNt 

T« = 

keNo 



yti = 



The corollary is of utility in evaluating det Fg 2 and det Fg 2 when 
det Fg[ and det Fg, are known, since 



det F 32 



and 



detF 



Sa 



yti = 

keNi 



7h = 
keNi 

are often considerably easier to evaluate than det Fg 2 and det Fg 2 . 
Of course similar remarks apply to Theorem VI. Frequently Theorem V 



374 THE BELL SYSTEM TECHNICAL JOURNAL, MARCH 1903 

also is useful in this respect. These results are generalizations of a theorem 
due to Bode. 

The following factorization theorem can be obtained by repeated 
applications of Theorem IV. 

Theorem VII: Let { (Ri , (R 2 , • • • , &„} denote any collection of disjoint 
subsets of (Ro such that U" =1 (R,- = (Ro . Then 



detF v = detFi.folldefFft.at 

i=2 
n 

det F 5 = del Fg.m, II del fW- 



7w = 

(k,i)e(Ri 

7ki = 

(k,i)e(Ri 



where <SU' = U;=l (R,- . 

Observe that det F 5 is expressed as a product of ordinary loop-differ- 
ences 4,5 when each <R,0' = 1,2, • • • , n) contains a single element. 

IX. SOME SIMPLE APPLICATIONS OF THE THEORY 

The first three examples relate to a specific vacuum-tube circuit that 
has been considered by Truxal. 6 He presents a detailed classical flow- 
graph analysis. 

9.1 An Application of Theorem I 

To illustrate the application of Theorem I we shall compute the 
driving-point impedance at port 2 of the circuit shown in Fig. 10 when 
port 1 is short-circuited. The pertinent linear incremental model is 
shown in Fig. 10. It is assumed that mi = M2 = 20, r pi = r p2 = 10, and 
Rli = RliR(Rli + R)~* = 200. 

Here we choose as the set ^ the two equations 

e a = Mi^i 

eb = M2002 

Let W = [e a , e b }' and «I> = [e i , e„o]', where the superscript t indicates 
transposition. Hence X = diag \jn , ps\- It is a trivial matter to show that 

440# fc + 4200 



w(0) = 



4.2tf* + 441 



Recall that C is defined here by 4» = C*F where the elements of W are 
treated as independent variables and port 2 is open-circuited. Similarly 
* = C*** when port 2 is short-circuited. An elementary analysis yields: 
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Fig. 10 — Circuit diagram for the example. 



c = 



-2.1/2, 



-2.1ft* 



C = 



4.2/?, + 441 




-7^(2/3 + 2.1) - 


420/3 


4.2/?* + 441 




-R k 




22/?* + 210 




-/?*(20/3 + 1) - 


200/i 



4.222* + 441 

gt(2g - 2.1) 
4.2/2* + 441 

-21/?* 
22/?* + 210 

**(g'-0 

22/?* + 210 J 



(23) 



22/? fc + 210 

The determinants det F 5 = det [1 2 - CX] and det F 5 = det [1 2 - CX] 
can be evaluated in a particularly simple manner by exploiting the fact 
that they are known at the outset to be linear fractional functions in 
R k . We find that 

12 - 840/?)/?* + 441 



det F 5 = 



det % = 



4.2/?* + 441 
(462 - 8400/3)/?* + 210 



(24) 



when mi = M2 = 20. Therefore 

u>(diag [20,20]) = w(0) 

= 200 
9.2 An Application of Corollary I 



22R k + 210 

det F 5 
det F 3 

(46.2 - 840/3)/?* + 21 
(88.2 - 840/3)/?* + 441 ' 



In this section let w(X) be the (port 1 to port 2) open-circuit voltage 
transfer function for the circuit in Fig. 10, and let $, »F, and <I> be 
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chosen as in the last section. Our objective here is to evaluate the trans- 
fer function in accordance with Corollary I which is obviously relevant. 
Since C and det F 5 have already been evaluated, consider the de- 
termination of the pertinent flow matrices A and B. By inspection of 
the circuit, A = [1,0]' and a simple calculation yields 

-2(R k + 210)1 



B 



2ft* 



4.2ft t + 441 



4.2ft* 



441 



It follows that 

det F 3 = det [1 2 - (C - AB)X] 
86.2ft fe + 441 



= det 



4.2ft fc + 441 

ft*. (40ft + 42) + 8400/3 
4.2ft fc + 441 



2R k - 20(420) 
4.2R k + 441 

(46.2 - 40p)R k + 441 



4.2ft* + 441 



(928.2 - 840/3 )R k + 441 + 160,000/3 



4.2ft* + 441 
Therefore from Corollary I, (24), and (25) 

840ft fc + 160,0000 



(25) 



w(diag [20,20]) = 



(88.2 - 840/3)ft* + 441 ' 



9.3 An Application of Corollary II 

Consider the model in Fig. 1 1 . Here let % denote the set of equations 

e a = Mi^i 

% = M26„2 

e k = R k i k 



0-/3)R 



/5R 



ik 



s 9= 



rr: 



e g , = -eK 



Fig. 11 — Linear model for the network in Fig. 10 when port 1 is short-circuited. 
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and let %i and $2 respectively denote the subsets 
e a = Mie ff i 

, e k = Rich- 

eb = V-&01 

Truxal shows 6 that det F S2 = 1 - [-0.2 + (40/21 )0\R k when the 
network parameters have the values given in Section 9.1. We wish 
to determine the corresponding expression for det Fg, using Corollary IT. 
By inspection of Fig. 11 

det F 5l p=° = 1, det F 52 ||— ° =l+±R k . 

Thus 

, . _ (88.2 - 840,3)fl* + 441 
det Fsi = *2R k + 441 

which of course is identical to the right-hand side of (24). 

9.4 A Flow Graph Demonstration of Corollary II 

Corollary II can be demonstrated by considering the flow graph in 
Fig. 12. Let the determinants of the loop-difference matrices for the 
branches Yi and Y 2 respectively be denoted by det F 5l and det F 5a . 
Straightforward evaluation shows that 

detF , = det [l m - Yi[d + K(l„ - Y 2 C 2 ) _1 Y 2 L]} 

detF S2 = det {1„ - Y 2 [C 2 + L(l m - YtCO^YiK]} 




Fig. 12 — Flow-graph relating to the validity of Corollary II. 
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where n and m respectively are the number of rows of Y 2 and the number 
of rows of Yi . 

The corollary implies that 

det \l m - Yjgi + K(l n - YATYJJj = det[l m - Y.CJ 
det { l n - Y 2 [C 2 + L( 1,„ - YA^YxK] } det [l n - Y 2 C 2 ] 

which can readily be verified independently with the aid of Lemma III. 

9.5 Foster's Results Via Theorem I 

As a final illustration of the theory we shall outline briefly an alterna- 
tive proof of Foster's well-known results concerning the readability 
of two-element-kind one-ports. 

Let an arbitrary passive RC one-port containing only resistors and q 
parallel-RC combinations with capacitance and resistance values re- 
spectively c k and {ac k )~ x (k = 1, 2, • • ■ , q) be characterized by the 
driving-point admittance function y(s + a), where a is real and posi- 
tive. Let g denote the set of capacitor volt-ampere constraints: i k = 
sc k e k (k = 1, 2, •■■,(() and consider the evaluation of y(s + a) in 
accordance with Theorem I (some of the flow matrices in Fig. 2 may 
not exist if a = 0). 

From Theorem I and Lemma I, y(s + a) can be expressed as 

y (s + « } = W(0) det[l g -,C 1 CC 1 ] (26) 

where d = diag [(ci)\ (c 2 )*, •■■ , (c,)*]. Observe that here -C and 
— C are passive open-circuit resistance matrices. A moment's reflection 
shows that A = B'. Hence — C = — C + D where D is a nonnegative 
definite matrix of unit rank. With the aid of (26) and the following 
theorem it becomes a simple matter to show that 

y(s) = go + «c» + £ - , *, 

fc=l S -+- D k 

where g Q , c M ^ 0; a k , b k ^ (k = 1, 2, ■ ■ ■ , q'). 

Theorem- VIII: Let P and Q denote two nonnegative definite hermitian 
matrices of order n; let Q have unit rank. Then 

P\ ^ n ^ 7> 2 ^ r 2 • • • ^ Pn £ r n 

where \p\ , p 2 , • • • , p n ] is the set of eigenvalues of P and {n ,r 2 , • ■ ■ , r„ J 
is the set of eigenvalues of (P + Q). Further, ifpj(j = 1,2, •■ -,n) is an 
eigenvalue of P of multiplicity m,-(j»y ^ 1), //ien pj is an eigenvalue of 
(P + Q) of multiplicity at least (my — 1) and at most (my + 1). 
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The discussion hi Appendix B shows that Theorem VIII is representa- 
tive of more general results that can be proved by remarkably simple 
arguments centering about Lemma IV. Eigenvalue inequalities of the 
type discussed in the Appendix are ordinarily deduced from the extremal 
properties of eigenvalues. 

APPENDIX A. 

Proof of Lemma I 

Since X X X 2 = X, Lemma I is a direct consequence of the following 
result. 

Lemma A: If D and E respectively are n X m and m X n ma- 
trices, det [l n + DE] = del [1 M + ED]. 

We prove first that the lemma is true when D and E are square matrices. 
Let D and E be p X p matrices. Then, if D is nonsingular, 

det [l p + DE] = det [D -, (l„ + DE)D] = det [1, + ED]. (27) 
If D is singular, it has a zero characteristic root, and hence there exists 
a positive number o- n such that D + al p is nonsingular for all real a 
satisfying < \ a \ < a a . Thus when < \ a \ < <tq , 

det [l p + (D + <rl p )E] = det [l p + E(D + *1 P )]. (28) 

Both sides of (28) are polynomials in a of at most degree p. Further- 
more these polynomials must be identical since they agree throughout 
the real interval (0, a ). Therefore (28) is valid when a = 0. 

Consider now the cases in which D and E are not square. Let 

p = m + n, 



D = 



m. 


ii 






n 


in 


[D 


°1 


n 
> 


E = 


[~E 


°1 





0_ 


m 




_0 


0_ 



m 



and let the symbol + denote a direct sum of matrices. Observe that 
det [l p + DE] = det [(1„ + DE) -j- 1,„] = det [1„ + DE], and that 
det [l p + ED] = det [(1,„ + ED) + 1„] = det [1,„ + ED]. 
This proves the lemma. 



APPENDIX B. 



On the Eigenvalues of a Sum of Matrices: an Application of Lemma IV 
In the following discussion M* and p(M) respectively denote the 
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complex-conjugate transpose and the rank of an arbitrary matrix M. 
Our principal result is 

Theorem A: Let \a i , a 2 , • • • , a„| denote the set of eigenvalues of the her mi- 
Man matrix A, where the first p eigenvalues vanish if A is of nullity p, and 

-(a p+1 )- x ^ -(a^y 1 ^ ^ -(a n y\ 

Let B denote a nonnegative definite hermitian matrix of unit rank and 
order n; and let 

. _ ,. det[U + s(A + B)] 
fc ~ S del [1. + sA] • 

Then — °o < k ^ °o and Me se£ 0/ eigenvalues of (A + B) can 6e written 
as {ci , c 2 , • ■ • , c„] where 

1. if < k < *> , the first p eigenvalues vanish and —(a p+i )~ ^ 
-(c p+1 )- J ^ -(a**) -1 ^ -(cp+2)" 1 ^ ••• ^ -(an) -1 ^ -(c)- 1 

2. if — co < & < 0, Me first p eigenvalues vanish and — (c p+ i) ^ 
-(ap+x)" 1 ^ -(c p+2 )- ! *^ -(flp+2)" 1 ^ ••• ^ -(*)"* ^ -(a^- 1 



3. if fc = 0, Me ,/rrsf (p + 1) eigenvalues vanish and — (a p+i ) ^ 
L (c p+2 )- 1 ^ -(a p+2 )- J ^ -(cp+s)" 1 ^ ^ -(c)- 1 ^ -(a„)" ] 

4. if k = °o ? Me ^rsf (p — 1) eigenvalues vanish and —(c p )~ ^ 
-(ap+0" 1 ^ -(cp+x)- 1 ^ -(a p+2 )-' ^ ••• ^ -(a,,)" 1 ^ -(c) -1 . 

Proof: 

The statements relating to the number of vanishing eigenvalues of 
(A + B) are a direct consequence of (i) p(A) = p(A + B) 
when < I fc I < 00, (U) p(A) = 1 + p(A + B) when k = 0, and 
(it*) p(A) = p(A + B) — 1 when k = ». 

Consider now the real rational function in ,s: 

r , = det [l n + s(A + B)] 
V{S) det[l„ + sA] 

It is well-known that there exist two unitary matrices P and Q such 
that PAP -1 = diag [ai , a 2 , • ■ • , a,,} and Q(A + B)Q _1 = diag [c x , c 2 , 
• • • , c n ] where the a,- and c, are real. Hence 

V(s) = II (1 +«;,-)( 1 + sa t y\ 

Consider the evaluation of y(s) in accordance with Lemma IV. Here 
B can be written as DD* where D is an n-vector. Thus 
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i—l-i 



y(s) = 1 + sD*[l„ + sA]"D 

= 1 + sD*P _1 [diag(l + soi , 1 + so* , ■ • ■ , 1 + sa„)] _, PD 
- 1 + sH*diag[(l + sai) _1 ,(l + so*)' 1 , •••,(!+ sa„) _1 ]H 

= l+±^il (29) 

,=i 1 + sen 

where H = PD = [/*i , /*■> , ■ ■ ■ , &»]'. The theorem is obviously true 
when ]C>"=i | hi | 2 = 0. Assume therefore that X)"=i I ^» P > 0- 

Observe that y'(s) > for all real s such that | y(s)\ < », and that 
the poles of y(s) are simple. Thus the interval on the real-axis of the 
s-plane between two adjacent poles contains one and only one zero of 
y(s). Note that when lim,^.«, y(s) = k > 0, the right-most critical 
point of y(s) is a zero and that this critical point is a polo if — =o < 
k ^ 0. Similarly the left-most critical point is a zero if k — <x> or — « < 
k < 0; it is a pole if ^ k < <x>. The inequalities stated in the theorem 
follow directly from these observations. The equal signs make provision 
for possible coincident nonzero eigenvalues of A and (A + B). Note 
that (29) implies that — °° < k ^ qo. 

The following corollary of Theorem A appears to be useful in the 
study of linear dynamical systems. 

Corollary A: Let A and B denote two hermitian matrices of order n, with 
A nonnegative definite. Let {ai , a 2 , • • • , a„\ in which a x 2£ a 2 ^ • • • ^ a n 
denote the set of eigenvalues of A; and let { C\ , c\ , • ■ ■ , c„) in which l\ ^ 
c-i ^ • • • ^ c n denote the set of eigenvalues of (A + B). Then 

St-fi ^'d [1 + £ * £ n] 

c, ^ d i+a [1 ^ i ^ (n - a)] 

where a = §(? + §), = h(f — s) in which f and s respectively are the 
rank and signature of B. 

Proof: 

Consider Theorem A and assume that A is nonnegative definite. Then 
either k > or k = °o [i.e., either p(A + B) = p(A) or p(A + B) = 
p(A) + 1]. Thus in either case, since the c< and a,- are nonnegative here, 

oi ^ ci ^ as ^ Co • • • g a„ ^ c„ (30) 

We may express B as XXi B, - Za=«+i B, where the B, are rank 1 
nonnegative definite hermitian matrices. It is certainly true that c, does 
not exceed the corresponding eigenvalue of A + X^"=i B, and is not less 
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than the corresponding eigenvalue of A — 2<=«+i B» • Hence by an 
a-fold application of inequalities of the type (30) we obtain the upper 
bound on c, stated in the corollary. If B in Theorem A were a nonposi- 
tive definite matrix the inequalities in (30) would be reversed, and 
hence a /3-fold application of the. inequalities suffices to establish the 
lower bound on c, . 

Our final result relates to eigenvalue multiplicities. 

Theorem B: Let A and, B be matrices of order n. Let A be similar to a di- 
agonal matrix and let the rank of B be unity. Then if a is an eigenvalue of 
A of multiplicity m(m ^ 1), a is an eigenvalue of (A + B) of multiplicity 
at least (m — 1). Further, if in addition A and B are hermitian, a is an 
eigenvalue of (A + B) of multiplicity at most (m + 1). 

Proof: 
Consider 

_ (s) = det[l n + s(A + B)] 
det[l„ + sA] 

Since B is of unit rank and A is similar to a diagonal matrix, the lemma 
can be used, as in the proof of Theorem A, to show that y(s) has only 
simple poles. This proves the first part of the theorem which is essen- 
tially equivalent to the statement that if the dimensionality of the null- 
space of [dl n — A] is m, then the dimensionality of the null-space of 
[dl n — (A + B)] is not less than (m — 1) [assuming that the Jordan 
form of A is diagonal and B is of unit rank]. It is not difficult to produce 
proofs of the result which are based on this interpretation. 

If A and B are hermitian, y'(s) is either positive for all real s such 
that | y(s) | < co or is negative for all such s (we assume that y(s) ^ 1 ), 
since B must be either nonnegative definite or nonpositive definite. Thus 
y(s) can have only simple zeros, from which the second part of the 
theorem follows at once. 
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